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CHAPTER I 


INTRODUCTION 

Multi-engine aircraft, with engines laterally separated from the 
aircraft longitudinal axis , experience a yawing moment whenever the en- 
gines are not producing equal amounts of thrust. This is nor- 
mally a liability, the severity of which depends on the: (1) Amount of 

thrust asymmetry involved; (2) Distance the engines are located from the 
aircraft centerline; (3) Direction of engine rotation for propellor 
driven aircraft; (4) Phase (s) of flight involved; (5) Aircraft config- 
uration, including asymmetrical fuel, weapons, or cargo loads as well 
as the position of flaps, slats, landing gear, etc.; (6) Other aircraft 
characteristics, including wing loading, total thrust available , center 
of gravity location, control effectiveness, stccompanylng loss of systems 
such as generators, hydraulic pumps, fuel pumps, or other systems which 
degrade the capability to properly assess and deal with the situation; 
(7) Environment (terrain obstruction, density altitude, field length, 
and weather); and, finally, (8) Actions of the pilot. 

While the number of variables involved is large, the options avail- 
able for corrective action are much more limited, and depend primarily 
on the phase of flight. Phases of flight can be classified as (Ref. 1) : 
(1) Nonterminal flight phases, and (2) Terminal flight phases. Non- 
terminal flight phases generally involve operations at relatively high 
airspeed, altitude, and Lift/Drag (L/D) ratios, and low angle of attack. 
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a, and angle of sideslip, 3. For the case of a significant thrust asym- 
metry (due typically to failure of one or more engines on one side of 
the aircraft), nonterminal flight phases such as cruise and loiter re- 
quire corrective actions that are much less urgent than for thrust asym- 
metry during terminal flight phases. The pilot can even temporarily re- 
duce thrust on the remaining englne(s) to eliminate the thrust asymmetry 
without Immediate fear of striking the earth or exceeding critical val- 
ues of a or B. There is time to attempt restarts or take other helpful 
actions without the immediate prospect of losing control. 

Terminal flight phases include takeoff, catapult takeoff, approach, 
wave-of f /go-ar ound , and landing. An engine loss can be critical during 
these phases, where the airspeed, L/D ratio, and altitude are relatively 
low, and a is relatively high. The corrective actions are urgent, and 
consist of (in sequence): 

(1) Selecting full power on all throttles. 

(2) Identifying the dead englne(s) and applying appropriate flight 
controls: 

(A) Rudder as necessary to maintain directional control. 

(B) Pitch control as necessary to maintain airspeed 
above the minimum for directional control (V^c) 
which is defined as the speed below which the eng- 
ine out yawing moment can no longer be controlled 
(some sideslip will exist) using the maximum avail- 
able rudder deflection. The pitch control must also 
keep a below stall, and the lift coefficient (C^) 
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high enough to prevent the aircraft from losing ex- 
cessive altitude and striking the earth’s surface. 
Note that the last requirement usually provides a 
serious conflict with the first two. 

(C) Aileron as necessary to counter rolling induced by 
the yawing moment, and to establish bank angles nec- 
essary to maintain directional control. 

(3) Making configuration adjustments to reduce the asymmetry, re- 
duce drag, reduce gross weight, reduce stores asymmetry, and 
so on. This might Include feathering a dead propeller, jet- 
tisoning external fuel tanks, changing flap settings, retract- 
ing the landing gear, closing cowl flaps or speed brakes, or a 
number of other possibilities depending on the aircraft. 

One of the primary reasons why thrust asymmetry can cause a crit- 
ical situation during terminal flight phases is because the dynamic 
pressure, which varies as the square of the airspeed, is relatively low. 
This results in reduced aerodynamic control effectiveness to counter any 
yawing moment produced by thrust asymmetry, since the yawing moment is 
relatively independent of dynamic pressure and does not decrease as air- 
speed is decreased. 

The use of vectored thrust to augment or replace aerodynamic con- 
trols is not a new concept. The German V-2 rocket of World War II used 
vanes located in the rocket exhaust to augment the aerodynamic control 
surfaces (Ref. 2). More recently, the British Hawker-Siddeley "Harrier" 
V/STOL jet uses thrust vectoring, and was placed in production in 1967. 
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Since then, a tremendous amount of material has been published concern- 
ing the use of vectored thrust in aircraft and missiles. Although most 
of the work In aircraft thrust vectoring has dealt with vectoring in the 
pitch plane for maneuver enhancement or improved V/STOL capability, the 
use of lateral vectoring has not been Ignored (Refs. 3-5). 

The purpose of this report is to investigate the use of laterally 
vectored thrust to counter the adverse effects of thrust asymmetry in a 
twin engine tactical jet, for the case of single engine failure during 
terminal flight phases. 

For this report , the linearized, small perturbation equations of 
motion were used to predict the steady-state control deflections re- 
quired for a single engine failure. Then, as explained in Appendices B 
and C, a full nonlinear six degree-of-freedom simulation was built on a 
hybrid computer to study transient responses. A hybrid computer was 
used for three main reasons: (1) The analog portion of the hybrid com- 

puter permits simulation of the real physical system as it actually per- 
forms (as a continuous system rather than discrete), and allows quick 
revisions to the system plant without extensive software modifications. 
The digital portion allows use of a modern digital control system to 
control the plant. The hybrid computer combines these capabilities and 
offers built-in analog-to-digital and digital-to-analog signal process- 
ing without having to interface analog and digital systems that were not 
specifically designed to interact. (2) It could easily Incorporate a 
cockpit, with analog displays and controls, at some future time, and 
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(3) It offered nearly unlimited availability and did not require use of 
a NASA Langley computer account. 

Finally, the simulator transient and steady-state responses were 
examined for different flight conditions and thrust vectoring param- 
eters, and the results compared to the analytical predictions. 

Unlike most American built, propellor driven twin engine aircraft, 
in which the left engine is the "critical engine," or worst case situat- 
ion for a single engine failure, jet powered tactical fighters do not 
have this limitation, and throughout this study only the case of a right 
engine out will be studied. All results will assumed to apply equal- 
ly to a left engine failure. 



CHAPTER II 


EQUATIONS OF MOTION 

The equations of motion for the aircraft can be derived from 
Newton's Second Law of motion. This is done in most aeronautical 
engineering textbooks; a thorough treatment is given in Reference 6. 

The usual assumptions used are; (1) the aircraft is a rigid body; 

(2) the mass of the aircraft remains constant for each particular 
flight condition examined, although the mass is different from one flight 
condition to another; (3) gyroscopic moments are not considered for 
any engines or rotating machinery on board the aircraft; and (4) the 
atmosphere is assumed fixed with respect to the earth. 

Under these assumptions the equations of motion can be written 
in the body axes as: 


Z - m(u + 
Z AF^ = m(v + 
Z AF^ = m(w + 
r = pi^ - 
E = ql^ + 
E AM^ = - 
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(2-1) 

( 2 - 2 ) 

(2-3) 

(2-4) 

(2-5) 

( 2 - 6 ) 
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The three kinematic equations involved are; 


ip = (q sin (j) + r cos (j))/cos 0 (2-7) 

0 = q cos (j) - r sin (f) (2-8) 

(J) = p + q sin 4> tan 0 + r cos (J) tan 0 (2-9) 


Euler angles ijj, 0, and (|) are defined in Figure 1; angles a 
and 3 in Figure 2. 

The sense of positive control deflections is also given by Figure 2. 
Note that the convention differs from NASA standard with regard to 
aileron deflection. 

The equations of motion were incorporated in the simulation through 
two steps. First, angular rate equations were rewritten from 
Equations (2-4), (2-5), and (2-6) as; 
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XZ 
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XZ 
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(r + pq) + Y“ 

X 



. ^Z 
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(2-10) 

( 2 - 11 ) 

(2-12) 


where body axis moments 


«x> 


My* and are obtained from the equations; 
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My = (2-14) 

™ P 

= qSbC + M_ (2-15) 

Z n ^ 

p 

where MXp, %p> ^ 2 p are the propulsive moments. 

Secondly, the force equations were transformed from body axes to 
wind axes, for two reasons: (1) Cl and Cp were available, not 
or C^; and (2) It would be difficult to program on the analog comput- 
er the Inverse trigonometric functions and square roots needed to solve 
for a and P when using velocity components In the body axes. The force 
equations In the wind axes are developed In Appendix A and can be writ- 
ten as: 


w 


m 


3 = p sin a - r cos a + 


w 


mV 


a = -(p cos a + r sin a) tan 3 + q + 


w 


mV cos 3 


(2-16) 

(2-17) 

(2-18) 


where F , F and F are defined in Appendix A. 

A. I 

w w w 

The equations of motion were completed by the expressions for ob- 
taining the velocity components in the body axes: 


u = V cos a cos 3 


(2-19) 


v = V sin 3 


( 2 - 20 ) 


w = V sin a cos 3 


( 2 - 21 ) 
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These components were also transformed into earth axes components 
by use of Euler angle transformations and then numerically integrated 
with the digital computer to obtain distances traveled over the earth, 
and altitude changes. All the necessary state variables for the simula- 
tion were then available, and were output to the analog strip chart 
recorders along with the control variables computed digitally by the 
Automatic Flight Control System (AFCS) subroutine. 

The numerical values needed to generate a particular solution are 
dependent on the flight condition involved. Aircraft physical charac- 
teristics and aerodynamic stability and control derivatives obtained 
from NASA are given in Appendix B. These data were used to define four 
straight and level, steady state flight conditions for dynamic analysis. 
The flight conditions are summarized in Table 1. The rational for 
selecting these flight conditions is given below. Although the follow- 
ing discussions deal with non-steady, untrimmed operation, the veloci- 
ties, angles of attack, and weight characteristics of each were used to 
specify the trimmed, straight and level cases given in Table 1. 

The conditions can be described as follows; 

CONDITION 1 - This is a "middle-of-the-road" first case. From the 
NASA data, the stability and control derivatives were approximated as 
linear functions in the a range from 0 to 20 degrees, so a = 10 
degrees and a reasonably low value of promised to be the most 

acceptable place to build and check the simulation. In addition, the 
available moments of inertia were explicitly given only for this weight 
distribution of fuel load (See Appendix B) . 
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CONDITION 2 - The weight was increased to include full Internal 
fuel, four radar missiles, full gun ammunition, and one externally 
mounted centerline tank, minus fuel necessary for engine start and run- 
up, 15 minutes of deck operations or taxi time, and an acceleration to 
79 m/s (154 KTAS). This speed is eight knots above the recommended 
field takeoff speed for a midrange center of gravity (26 percent MAC). 

It also approximates a catapult takeoff "end airspeed," which is the 
airspeed attained by the aircraft at the end of the catapult power 
stroke. Note that the thrust value given in Table 1 is for steady-state 

flight and is not representative of takeoff thrust. The inertias for 

this and the following conditions were calculated in accordance with the 
discussion of Appendix B. 

CONDITION 3 - For the same weight and moments of inertia as Con- 
dition 2, the Initial speed was set at five knots below the computed 

takeoff speed. This represents the loss of airspeed normally encounter- 

ed during a single engine failure while taking off at high gross 
weights, if corrective action is not prompt and correct. It also rep- 
resents a common operating point for such an aircraft when involved in 
the primary mission of air-to-air combat, where a single engine failure 
at relatively high a and low dynamic pressure can cause loss of direc- 
tional control and subsequent Spin entry. 

CONDITION 4 - For the same a as Condition 3, but at a reduced 
weight representative of a landing approach condition, Condition 4 sim- 
ulates a single engine failure during a wave-off /go-around maneuver. 



CHAPTER III 
STATIC CONTROL 

In order to analytically determine control deflection requirements 
for a single engine failure without resorting to analog or digital com- 
putation, the nonlinear Equations (2-2), (2-4), and (2-6) must be 
linearized. If this is done, it is possible to quickly obtain solutions 
which can be used to direct the research effort, or to verify the 
computer solutions obtained from the nonlinear equations. The 
linearized, steady-state equations for the lateral motion were used. 
These equations were obtained from the general equations by incorporat- 
ing two assumptions: (1) the aircraft motion is restricted to small 

perturbations from a reference condition of symmetric steady flight 
with no angular velocities; and (2) since the perturbations are assumed 
small, the products of perturbations can be neglected. 

The form of the linearized equations defining steady sideslip 
is: 


■Fy -mg cos 0 • (}) 

__P 

T[S 




(3-1) 
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6 

r 



(3-3) 


Since the thrust vector lies in the X-Y plane of the aircraft 

and was not vectored in the pitch plane, =0. When the thrust is 

P 

asymmetric and non*vectored = 0, because the non-vectored thrust 

line is assumed parallel to the X axis. For Equations (3-1) through 
(3-3), all the parameters except (|), 3, , and 6 are known for each 

IT 3 . 

trimmed initial flight condition. 

The APCS was designed to use bank angle and rudder and aileron 
deflection to achieve 3=0. Equations (3-1), (3-2), and (3-3) 
can be rewritten to represent this condition, which leaves three linearly 
independent equations in three unknowns ((J), 6 , 6 ) to be solved 

IT 3 

simultaneously. These equations can be written in a state-space equation 
of the form y = Ax, where 


E. 



^b 
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Then x ” ^ 

Calculations with an engine out for both the non-vectored case and the 
case when the remaining thrust is vectored 5.7 degrees (through the 
aircraft center of gravity) resulted in control deflections as shown 
in Table 2. 

The calculations indicate that vectoring the thrust will eliminate 
the requirement for rudder and aileron deflections, and also reduce the 
required bank angle slightly. Although the elimination of aileron 
deflection is hardly significant, since the greatest deflection required 
for the rK)n-vectored case is still less than 1 degree, the elimination of 
13.5 degrees of rudder deflection is a significant achievement 
(maximum rudder deflection possible is 30 degrees) . 



CHAPTER IV 


DYNAMIC CONTROL 

Once the basic aircraft simulator had been built and the fundamen- 
tal modes of longitudinal and lateral motion verified, it then became 
necessary to control the dynamic motion of the aircraft. If subjected 
to thrust asymmetry without an active controller, the aircraft flight 
path would quickly diverge from the initial trimmed state, and the 
aircraft would "crash." In order to control sideslip, bank angle, and 
all other parameters of interest, it was necessary to provide control, 
either by having a pilot manually supply control inputs through an an- 
alog cockpit connected to the simulator, or by programming an automatic 
controller to perform the same tasks. For this early stage of the sim- 
ulator's maturity, the automatic controller was chosen in order to by- 
pass the additional engineering a cockpit would have required, and in 
order to remove the variability of human response. A block diagram of 
the automatic controller is given as Figure 3. The automatic control- 
ler is labeled "AFCS" (Automatic Flight Control System) , and incorpor- 
ates a Stability Augmentation System (SAS) with inner loops to augment 
basic aircraft stability and control, and an Autopilot with outer loops 
for flight path control. As shown, the SAS provides phugoid damping, 
Dutch roll damping, roll rate damping, and sideslip control, while the 
Autopilot provides altitude hold, heading control, and wing leveling. 

As explained in Appendix C, the Pacer 100 Digital Processor was 
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used to accomplish the SAS and Autopilot functions, as well as to pro- 
vide the additional functions of engine dynamics control and Thrust Vec- 
tor Control (TVC) shown In the "Thrust modes" block of Figure 3. 

Since the TVC mode Is of primary Importance for this report, an In- 
troductory discussion Is In order. 

If thrust Is to be vectored, there are two fundamentally different 
types of control that can be used. The first Is a fully active, closed- 
loop controller that can be used to provide or augment directional stab- 
ility, as outlined In Reference 4. That report employed a non-specific 
two-dimensional TVC nozzle, which Is a typical approach for a fully ac- 
tive system. The second type Is an open-loop controller that Is not 
used unless activated by the pilot or another control system In response 
to an engine out. An open- loop controller can make use of a simpler 
vane system that Is external to the exhaust jet, and which does not In- 
terfere with the exhaust jet during normal twin engine operations. 
Research Is currently being conducted In this area (Ref. 7), and one 
possible scheme employs two vanes deflected simultaneously and equally 
when activated (see Figure 4). In this way, no matter which engine 
failed, thrust from the remaining engine would be vectored through the 
aircraft eg. Also, In case of Inadvertent activation, the net result 
would only be a small percentage loss of longitudinal thrust component 
(0.5 percent for a thrust vector angle of 5.7 degrees, not Including any 
efficiency losses). This study used an open- loop controller to drive an 
external vane system, because of the simplicity and inherent safety. 

The Thrust Vector Controller was incorporated in the simulation 
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as shown In Figure 3. Three decisions were made for each digital 

5^ going from the AFCS to the "Thrust mode" controller. These 

command 

were open-loop decisions; i.e., independent of AFCS inputs or outputs 
or of the parameters generated by the aircraft simulation. These de- 
cisions were made by the person operating the simulation, through three 
switches located on the Digital Processor console, and shown schematic- 
ally in Figure 3. The experimenter could select various combinations 
of the thrust modes shown. One switch, when activated, set the right 
engine thrust equal to zero. A second switch either activated TVC logic, 
or else bypassed it in order to keep the thrust conventionally non- 
vectored. A third switch controlled selection of the thrust "command 
double," which will be explained later. If the experimenter activated 
TVC logic, the computer required three additional pieces of information; 
(1) The desired reaction delay time, from the moment of single engine 
failure until the vanes begin to deflect; (2) the actuation interval, or 
time that it takes the vanes to fully deflect the thrust vector once the 
vanes begin to move; and (3) the desired thrust vector angle, 6^, which 
will vector the thrust through the aircraft eg. This angle can be cal- 
culated from the aircraft geometry as shown in Figure 4. 

This report does not address the efficiency of the vanes in turn- 
ing the thrust vector. Based on results reported in Reference 7, it is 
assumed that a thrust vector angle of 5.7 degrees for this aircraft can 
be achieved with vane deflection angles less than 10 degrees. 

At this point, since a known constant 6^ is assumed, the parameters 
involved in investigating the engine out case are five-fold: 
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(1) Which flight condition Introduced in Chapter III is being consider- 
ed; 

(2) Whether or not the remaining thrust will be vectored in response to 
an engine out; 

(3) If vectored, the reaction delay time used, as well as 

(4) The actuation Interval used, and 

(5) Whether or not the thrust command double was used. 

After obtaining baseline data runs for each of the four flight con- 
ditions without thrust vectoring, the TVC was used on all remaining 
runs. Reaction delay times used were two and four seconds, since only 
one second is an unrealistically short time for a pilot to identify an 
engine out condition and activate a control system. Actuation Interval 
times used were one and two seconds, since four or nwre seconds to move 
a vane 10 degrees or less is uncharacteristically slow. Each combina- 
tion mentioned above was run both with and without the thrust command 
double. 

The reason for the thrust command double is as follows; If only 
the altitude error and phugold damping requirements are used to generate 
a thrust command, the thrust response is too slow, and excessive alti- 
tude is lost when an engine is cut. Therefore, a mode was added which, 
when selected, doubled the commanded thrust of the remaining engine as 
soon as the reaction delay time expired. This commanded thrust was 
still subject to the engine dynamic limitations, but it quickly increas- 
ed thrust from the single engine in order to equal the initial value of 
total thrust from both engines. This type of response is closely akin 
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to a normal pilot's reaction to an engine out case, when all throttles 
are Immediately advanced to full power, before Identifying which engine 
failed. 

For the five parameters mentioned, there are a total of four 
baseline non-vectored cases to consider, and 32 vectored cases, for a 
total of 36. These test cases are denoted by an ’*X" In Table 3, the 
test matrix. 

The results will be presented In Chapter V by first examining the 
baseline non-vectored cases for each of the four flight conditions, and 
then Investigating the effect of parameter variations. Use of the com- 
mand double Is Investigated first, followed by a systematic variation of 
actuation Interval and reaction delay times. 



CHAPTER V 


RESULTS AND DISCUSSION 

Figures 7 through 14 show the time histories for the baseline non- 
vectored cases of the four chosen flight conditions. The transient peak 
values and steady-state values are summarized in Table 4(a), and the 
steady-state values of <|), 6^, and 6^ agree well with the calculated 
values in Table 2. Conditions 3 and 4 have the largest peak and steady- 
state values, and, except for altitude perturbations (Az), the steady- 
state values of interest for Conditions 3 and 4 are equal, and peak val- 
ues are very similar. This was also true for the thrust vectored cases. 
Steady-state values of <|), 3, 6^, and 6^ for the thrust vectored 
cases are shown in Table 4(b), and also agree well with the calculated 
values in Table 2. 

The only figures for vectored thrust that are included in this re- 
port are for Condition 3, in order to avoid duplication of Condition 4, 
and in order to represent the worst overall condition to have to con- 
trol. Also, Condition 3 figures are restricted to the case for the 
shortest vectored response time (2 sec reaction delay/1 sec actuation 
interval) and the case for the longest vectored response time (4 sec 
reaction delay/2 sec actuation interval). Intermediate cases were re- 
dundant . 

Since the peak transient values of lateral variables for Condition 
4 are slightly greater than for Condition 3, with steady-state values 
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being equal, the problem of providing lateral mode stability and control 
is just as great at the reduced weight and speed representing a landing 
approach (Condition A) as for the higher gross weights associated with a 
takeoff (Condition 3). Of course, the descent rates and altitude losses 
are not as great for the reduced weight condition, but any given propul- 
sive yawing moment can produce higher values of p and r in the con- 
dition with the lower weight and moments of inertia. 

As explained on page 17 of Chapter IV, it was necessary to use the 
thrust command double to prevent unacceptable altitude losses in all the 
flight conditions. Examination of Condition 3 (which has the highest 
combination of weight and a) emphasizes this necessity: Figure 11 shows 

an altitude perturbation of 450 m (1476 ft) for the baseline non- 
vectored case of Condition 3, vdiich does not use the thrust command 
double. In fact, the maximum available thrust of the remaining single 
engine at military power [67,500 N (15,175 lb)] is just adequate to can- 
cel the descent rate incurred from the single engine failure, and the 
simulated aircraft maintains level flight 400 m (1312 ft) below the 
initial altitude, because there is no excess thrust to regain the alti- 
tude lost. When the thrust command double is used with Condition 3, as 
shown in Figure 17, the altitude perturbation is only 50 m (164 ft). 

The quicker engine response is important in providing safety margins and 
pilot acceptance, and should be incorporated in any investigation of the 
longitudinal or lateral variables. Use of the command double caused 
greater perturbations in the lateral variables, which can be seen by ex- 
amining Figures 16, 18, 20, and 22 (the data are summarized in Table 5). 
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This increase averaged 10 percent for the worst case of the study, but 
this is acceptable in light of the altitude considerations. 

Although use of thrust vector control can eliminate the need for 
lateral steady-state aerodynamic control deflections, it is important to 
also consider the effects of TVC on longitudinal performance. As men- 
tioned on page 14 of Chapter IV, 5.7 degrees of vectoring does decrease 
the longitudinal thrust component by 0.5 percent. Comparison of Figures 
11, 15, and 19 (for which the appropriate variables are summarized in 
Table 6) shows that this decrease is more than offset by benefits de- 
rived from thrust vectoring: (a) Eliminating the lateral steady-state 

aerodynamic control deflections reduces the aircraft drag by approx- 
imately 1.4 percent; and (b) The maximum bank angle is reduced from 30 
degrees to 14.5 degrees, and comparison of Figures 12 and 16 shows that 
the area under the curve for bank angle is approximately one-fifth as 
great as for the non-vectored case. This translates into more effective 
lift and less altitude lost as a result of bank excursions. 

Therefore, the net result of TVC on longitudinal variables for this 
aircraft is to enhance performance. An equivalent 1.4 percent loss in 
longitudinal thrust component occurs for a thrust vector angle of 9.6 
degrees, implying that vector angles less than that offer potential for 
improving longitudinal performance. 

Finally, the effect of larger reaction delay and actuation interval 
times is to increase the peak values of the lateral variables while de- 
grading longitudinal performance. This statement can be generalized to 
Include all four flight conditions. Of note here is that the 6^ for 



22 


Conditions 3 and 4 reaches the physical limit of 20 degrees deflection 
for both the non-vectored case, and the vectored case with a 4 second 
reaction delay. 6^ saturated for approximately 16 seconds in the 
non-vectored case and 3.5 seconds In the vectored case. By using a 
reaction delay of 2 seconds with the TVC, control surface saturation is 
totally avoided . 



CHAPTER VI 


CONCLUDING REMARKS 

The purpose of this report was to investigate the use of laterally 
vectored thrust to counter the adverse effects of thrust asymmetry in a 
twin engine tactical jet, for the case of single engine failure while 
operating at relatively low airspeeds and L/D ratios typically encoun- 
tered in landings, takeoffs, and air combat maneuvering. There were 
several results of the research effort: (1) The linearized, steady- 

state, lateral equations of motion (used to analytically determine the 
control deflections required for the engine out case) closely predicted 
the results provided by flight simulation; (2) The non-linear, six 
degree-of-f reedora simulator that was built on a hybrid computer accu- 
rately represents the modeled aircraft for values of a from 0 to 20 de- 
grees (based on comparisons of the simulator’s fundamental modes of lon- 
gitudinal and lateral motion with known values for the data used), and 
it will provide a flexible tool for further studies of vectoring the 
thrust laterally, vertically, or in combination. The simulator has the 
potential to include a piloted cockpit, and to operate at greatly in- 
creased values of a; (3) For the flight conditions and aircraft geom- 
etry Investigated, laterally vectoring the thrust 5.7 degrees, within 
three seconds of engine failure, required peak rudder deflections less 
than half that needed when thrust was not vectored, and reduced steady- 
state deflection from 13.5 degrees to zero. Thrust vectoring similarly 
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reduced maximum bank angles and heading errors by less than half, while 
preventing saturation In aileron deflection. 

From these results, several conclusions can be drawn for the case 
of single engine failure in the aircraft studied: 

(1) Lateral thrust vectoring can provide significantly decreased peak 
values of rudder and aileron deflection, heading error, bank angle, and 
yaw rate. It can eliminate the requirement for significant steady-state 
rudder deflection. 

(2) For relatively small thrust vector angles, the decrease in the lon- 
gitudinal thrust component caused by the change in the thrust direction 
can be offset by improved aerodynamic efficiency of the airplane because 
the aerodynamic controls return to their neutral positions. This as- 
sumes that the thrust vectoring device is perfectly efficient in chang- 
ing the thrust direction. 

(3) The reaction delay time (from the moment of single engine failure 
until the thrust begins to be vectored) is an important parameter which 
ideally is zero (no delay). If the reaction delay is too long, lateral 
aerodynamic control saturation may be encountered, even though the vec- 
toring will eventually reduce steady-state deflection requirements to 
zero if the aircraft does not go out of control. 

(4) Single engine failure during the relatively low weight landing ap- 
proach flight condition can present lateral stability and control prob- 
lems equal to or worse than those encountered during takeoff at higher 
weights. 

Based on the results and conclusions, laterally vectored thrust 
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appears to offer a number of additional potential advantages: (a) The 

ability to augment directional stability through a closed- loop, active 
controller, thereby decreasing vertical stabilizer size, which offers 
weight savings and reduced drag; (b) The ability to generate anti-spin 
moments, thereby lessening the potential for aircraft loss from out-of- 
control flight; (c) The ability to generate side-force components with- 
out bank angles, which could Improve an aircraft's capability to handle 
crosswinds during landing approaches to an airfield, and to make lineup 
corrections during approaches to an aircraft carrier; and (d) Mainte- 
nance of controllability and directional stability at significantly in- 
creased a and lower dynamic pressure could play a large role in safely 
lowering operating speeds during terminal flight phases or air combat 
maneuvering, either one of which can enhance the effectiveness of a tac- 


tical jet aircraft. 



APPENDIX A 


DEVELOPMENT OF FORCE EQUATIONS 
From Newton's Second Law of motion: 


I AF m 


dt 


In the wind axes: 


mg + F + F 
2w — a — p 
w *^w 
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■jr + (0 V V 
dt — w ^ — 


where from References: 


= pi + q_j+rk 
"w *^w- ~w 
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= p cos a cos 3 + (q ” a) sin 3 + r sin a 
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where 


(X 


= -(p cos a + r sin a) tan 3 + q + 


-w 


mV cos 3 


F =» mg cos 3(cos 0 cos ({> sin a - sin 0 cos a) 

X 

w 

+ mg sin 3 cos 0 sin 4* + ^S(Cy sin 3 - cos 3) 


+ cos 3 


cos a + ^2 
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APPENDIX B 


CHARACTERISTICS Of AN ADVANCED FIGHTER 

Aircraft physical characteristics and aerodynamic stability and 
control derivatives were obtained from NASA for a representative swept- 
wing, Mach 2 class twin engine fighter at one operating weight. The 
physical characteristics are: 

S = 56,48 m^ « 34,574 kg-m^ 

c * 4.86 ro 1^-225.900 

b « 13.10 m 

m - 16,280 kg Ij ^2 ▼ 

W « 159,652 N - 7.0 m (longitudinal distance 

from exhaust nozzle 4ixit 
plane to aircraft eg) 

Zy. * 0.7 m (lateral distance from 
center of exhaust nozzle to 
aircraft centerline) 

Basic aircraft weight without fuel is approximately 122,500 N 
(27,539 lb). 

Since moments of inertia are proportional to the mass involved, 
and proportional to the square of the radius arm involved, it was 
assumed that the relative distribution of any fuel loads different 
from that originally given would be the same as the original distribu- 
tion, though the mass would differ. Scaling the moments of inertia 
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to the relative masses chosen in Chapter II for the four different 
test conditons results In the moments of inertia for each test condi- 
tion as shown in Table 1. 

The configuration was clean, with landing gear, flaps, and speed 
brakes retracted. Since aerodynamic data for the case of landing 
gear and flaps extended was not explicitly available, the clean 
configuration was retained for the study rather than having to make 
aerodynamic assumptions to extend the available wind tunnel data. 

The stability and control derivatives do not include variations 
due to speed effects since the given data were for a specific Mach 
number, and the speed band involved is only 24.8 m/s (48.2 knots). 

The study is confined to the flight regime below stall a, where 
variations in derivatives are essentially linear. As a result, 
simplifying approximations were made by writing the derivatives as 
linear functions of a. The equations used for the coefficients were: 
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D D D- L Dp r Yo 

o 0 3 

L r 


C =Ca-fC 6+~Cq 
m m mp e 2V m 

a 6 q 

e ^ 


(B-1) 

(B-2) 

(B-3) 




(B-4) 



3 1 


C„ =C„3+C« 6„+C„ + -;jjrlC„ p + C 


•A ” " ■ ''i. "r a •2V 

3 6^ ^ ^ P 


C «C3 + C 6+C 6+^(Cp + C 

n iig *" tig a 2V I n 

r a 


(<=<ipP + 

( V * v) 


(B-5) 


(B-6) 


where (all values per radian): 


C„ « -0.0158 + 0.0615a 
a 


» 0.1 + 0.8881a 
r 


C- - 0.0115 - 0.0327a 
6 

r 


C. - 0.058 - 0.1047a 
a 


C = 0.165 - 0.525a 
"p 


C = -0.286a 
n 

P 


C = 0.00573 - 0.0274a 
a 


All other derivatives are constant, with the following values 


Cj^ = 3.466 


C = -10.6 

IQ 


(L » 0.544 
e 


s. 


-0.91 


Cj, - 0.337 

S, 


S, 


- 0.174 



= 0.088 


0. 129 




C = -0.347 



€ = -0.693 

0 

e 



C = -0.43 
n 

r 


The values of and for a = 0 are: 


C. = 0.025 

Li 



The control deflection limits are: 

6^ = -35, +15 degrees 
6^ = 20 degrees 
6^ = 30 degrees 
(S =5.7 degrees 


-0.272 


-0.084 


0.0175 


T = 67500 N per engine 



APPENDIX C 


DESCRIPTION OF EXPERIMENTAL FACILITIES 
The hybrid computefr used in this study belongs to the NASA Langley 
Research Center’s Aerospace Controls Research Laboratory. It is an 
Electronic Associates, Inc., (EAI) model, incorporating a Pacer 100 
Digital Processor with 32,000 sixteen-bit words of core memory, a 
681 Parallel Analog Processor, and a single fixed-head disc storage 
unit capable of storing 360,448 sixteen-bit words (Ref. 8). The 
Pacer 100 memory cycle time is 1.0 microsecond, with a subtract or 

divide execution time requiring 2.0 or 6.6 microseconds, respectively. 

/ 

Peripheral Input /Output (I/O) devices include: INPUT - Paper tape, 

cassette magnetic tape, cards, and interactive terminal; OUTPUT - Paper 
tape, cassette tape, line printer, interactive terminal, X-Y plotter, 
and eight-channel strip chart recorders. See Figure 5. The sixteen bit 
I/O Bus allows communication with the Pacer 100 at over 555,000 words per 
second. The system uses a Real-Time Clock Unit, initialized from the 
I/O Bus, to synchronize program operation between the Digital Processor 
and Parallel Analog Processor. The Parallel Analog Processor was used 
to program the nonlinear equations of motion for three of the six degrees 
of freedom (angular displacements) , and the Pacer 100 Digital Processor 
was used to calculate the remaining three degrees of freedom (linear 
displacements) as well as to set the analog potentiometers that acted as 
coefficients for both the calibration checkout runs and the real-time 
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data gathering runs. In addition, the Digital Processor was used to 
implement engine dynamic characteristics, thrust vectoring logic, and 
AFCS logic for inner-loop stability augmentation and outer-loop auto- 
pilot functions. The simulation proceeded through the steps outlined 
below and shown in the simplified flow chart of Figure 6. First, 
aircraft physical characteristics and aerodynamic stability and control 
derivatives (obtained from NASA for a representative Mach 2 twin engine 
fighter) were read into Digital Processor storage arrays from a punched 
paper tape. Test point parameters and initial conditions (ICs) , which 
could be varied, were also read into storage arrays at the same time. 
Then, the Digital Processor used the appropriate array values to set the 
variable potentiometers of the 681 Parallel Analog Processor. 

Calibration checks were then accomplished for a known operating trim 
point. If the calibration checks were passed, the main operating pro- 
gram was loaded into the Digital Processor from punched cards, which 
simultaneously activated the analog hardware and the repetitive digital 
loop. Once the Parallel Analog Processor was activated from the initial 
conditions, it operated on the aircraft equations of motion dealing 
with angular displacements and angular rates, at a speed of 12 million 
Equivalent Operations Per Second, which is significantly faster than the 
Digital Processor’s speed of 300,000 operations per second (based on 
the Gibson mix of fixed-point/floating point arithmetic and logic 
operations). Equivalent Operations Per Second represents the speed of 
a digital computer that would perform the same computations in the same 
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time and with the same accuracy (to three decimal places in this case) 
as the analog computer. 

At fixed time intervals dictated by the digital controlling pro- 
gram, selected variables from the analog side were sampled, converted 
to digital format, and input to the Digital Processor to be used in 
digital calculations. The digital Central Processor Unit (CPU) had to 
handle three jobs during program operations: (1) direct the overall 

program in terms of sampling analog channels at appropriate times, 
outputting digital signals to be converted to analog format at correct 
times, and updating aerodynamic coefficients; (2) process the linear 
displacement equations, which required analog inputs; and (3) process 
the control laws. Four basic tasks were accomplished in the Control 
Laws subroutine: modeling engine dynamics in terms of thrust limits 

and time responses; modeling Thrust Vector Control (TVC) logic in terms 
of delay times and rate of vectoring; SAS calculations; and Autopilot 
calculations. Analog inputs were needed for the last two tasks. 

Once the Control Laws subroutine was completed and appropriate 
control Inputs had been computed, these inputs were converted from 
digital format to analog where necessary, and then fed back to the 
digital or analog processor to form the closed loop system. Two eight- 
channel analog strip chart recorders were used to graph the control 
deflections and longitudinal and lateral variables as functions of 
real operating time. 

The digital CPU scheduled the three jobs mentioned previously by 
using a Priority Interrupt system that gave the Control Laws subroutine 
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the highest priority, and made the CPU available to process the control 
laws every 40 milliseconds. Second priority was given to processing 
the linear displacement equations, which occurred every 100 ms, and 
the remaining time was left for overall program execution. This 
prioritizing prevented a situation where the control laws would be 
left with too little time to complete necessary calculations, which 
could lead to system stability problems. With the control laws given 
first priority, any saturation occurring in that portion of the loop 
would stop the program and alert the experimenter to the problem. 

As functions were added to the Control Laws subroutine, an oscilloscope 
was used to check the time being taken for the subroutine calculations, 
to insure that limits were not being approached. In its final form, the 
Control Laws subroutine required approximately 10 milliseconds for 


completion. 
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TABLE 1. - CHARACTERISTICS AND PARAMETERS FOR EACH FLIGHT CONDITION TESTED 





Condition 1 


Condition 2 

Condition 3 

Condition 4 

Weight 

, W, N 

(lb) 

159 652(35 891) 

201 

950(45 400) 

201 950(45 400) 

150 350(33 800) 


2 

-m 


34 574 


43 734 

43 734 

32 560 




225 900 


285 750 

285 750 

212 738 




253 540 


320 713 

320 713 

238 768 

I „ ^ 



-13 558 


-17 150 

-17 150 

-12 768 

xz 








a^, deg 


10 


15 

17.5 

17.5 

Mass, m, kg 


16 280 


20 593 

20 593 

15 332 




.58 


.86 

1.00 

1.00 

S 



.13 


.27 

.35 

.35 

Thrust 

, T, N 

(lb) 

35 262(7 927) 

59 

880(13 462) 

67 484(15 171) 

50 241(11 295) 

Altitude, h. 

m 

sea level 


sea level 

sea level 

sea level 

V, m/s 

(knots) 

87.27(170) 


79.07(154) 

72.41(141) 

62.48(121.5) 

q, N/m 

2 


4 669 


3 833 

3 214 

2 393 



TABLE 2 - REQUIRED STEADY STATE CONTROL DEFLECTIONS FOR RIGHT ENGINE OUT AND /S = 0 

Condition 1 Condition 2 Condition 3 Condition 4 

Unvectored Vectored Unvectored Vectored Unvectored Vectored Unvectored Vectored 


Sideslip, 3, deg 

0 

0 

0 

0 

0 

0 

0 

0 

Bank, 4 ), deg 

-1.4 

-1.3 

-1.9 

1 

0 

-2.2 

-2.0 

“2.2 

-2.0 

deg 

4.8 

0 

10.1 

0 

13.5 

0 

13.5 

0 

6 , deg 

-0.7 

0 

-0.9 

0 

-0.7 

0 

• 

O 

1 

0 
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TABLE 3. - TEST MATRIX 

Condition 1 Condition 2 Condition 3 Condition 4 


Basic Nonrvectored 

x” 

x» 

x^ 

x» 

Vectored 

thrust 





Reaction 

Actuation 





delay 

interval 





sec 

sec 





2 

1 

X 

X 

x" 

X 

2 

1® 

X 

X 

x" 

X 

2 

2 

X 

X 

X 

X 

2 

2® 

X 

X 

X 

X 

4 

1 

X 

X 

X 

X 

4 

1® 

X 

X 

X 

X 

. 4 

2 

X 

X 

x» 

X 

4 

2® 

X 

X 

x'’ 

X 

^Includes 

thrust command double 




^Included 

in List of 

Figures 






TABLE 4 . - SIMULATION TKffi HISTORY VALUES 
(a) Basic non-vectored flight conditions 



Condition 1 

Condition 2 

Condition 3 

Condition 4 


Peak 

Steady-state 

Peak 

Steady-state 

Peak 

Steady-state 

Peak 

Steady-state 

V, m/s 


87.27 


79.07 


72.41 


62.48 

a» deg 

12 

10 

18 

15 

20 

17.5 

20 

17.5 

q, rad/s 

-.007 

0 

-.01 

0 

.03 

0 

.035 

0 

z(-6), m/s 

8.2 

0 

14 

0 

22 

0 

21 

0 

Az(-Ah), m 

220 

0 

350 

0 

450 

400 

375 

0 

T, N 

40 000 

35 262 

67 500 

59 880 

67 500 

67 484 

58 000 

50 241 

p, rad/s 

.03 

0 

.085 

0 

.102 

0 

.115 

0 

<[>, deg 

2.6 

-1.4 

11.5 

-1.89 

-30 

-2.18 

-30 

-2.18 

r, rad/s 

.02 

0 

.035 

0 

.075 

0 

.09 

0 

3 , deg 

-1.1 

0 

-2.7 

0 

-3.3 

0 

-3.0 

0 

deg 

5.4 

4.8 

11.1 

10.2 

14.4 

13.5 

16 

13.5 

6 , deg 

8L 

-4.4 

-0.6 

-16 

-.6 

-20 

-.6 

-20 

-.6 

deg 

4.0 

2.0 

17.5 

5.6 

48.6 

5.6 

64.2 

5.6 



<l>, deg 
3, deg 
<Sj., deg 
6 , deg 


TABLE 4 . - Concluded. 

(b) Vectored thrust steady-state values. 


Condition 1 
-1.29 
0 
0 
0 


Condition 2 
-1.72 
0 
.15 
0 


Condition 3 
-1.95 
0 
0 

.2 


Condition 4 
-1.95 
0 
0 

.2 


N> 



TABLE 5 . - EFFECT OF COMMAND DOUBLE ON VECTORED THRUST LATERAL VARIABLES 


Condition 3 


Condition 3 


Condition 3 Condition 3 


2 sec reaction 2 sec reaction 4 sec reaction 4 sec reaction 
1 sec actuation I sec actuation 2 sec actuation 2 sec actuation 
No coim&and double Command double No command double Command double 



Peak 

Steady-state 

Peak 

Steady-state 

Peak 

Steady-state 

Peak 

Steady-state 

p, rad/s 

-0.132 

0.0 

-0. 132 

0.0 

-0.132 

0.0 

-0. 132 

0.0 

deg 

14.5 

-2.9 

14.5 

-2.9 

21.7 

-2.9 

23.6 

-2.9 

r, rad/s 

.04 

0 

.04 

0 

.053 

0 

.06 

0 

3, deg 

-3.0 

0 

-3.0 

0 

-3.2 

0 

-3.2 

0 

6^, deg 

6.6 

0 

6.6 

0 

7.2 

0 

7.8 

0 

6 deg 

a 

-17.0 

.2 

-18.5 

.2 

-20 

.2 

-20 

.2 

, deg 

8.6 

2.9 

8.6 

2.9 

15.7 

2.9 

17.5 

2.9 



TABLE 6. 



Condition 3 
Non-vectored 


Variation in V> m/s 

64.5-80.0 

Maximum a, deg 

20 

Maximum S , deg 

-10.8 

Maximum z, m/s 

22 

Peak/S teady-st ate Az, m 

450/400 

Steady-state T, N 

67 484 


i 

I 

i 

i 

i 

I 


ON LONGITUDINAL VARIABLES 


Condition 3 

Condition 3 

2 sec reaction 
1 sec actuation 
No command double 

4 sec reaction 
2 sec actuation 
No command double 

64.0-73.5 

65.0-75.5 

20 

20 

-10.2 

-10.5 

15 

17.75 

370/350 

375/350 

67 484 

67 484 
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Figure 1. - Definition of Euler angles. 







(a) Right engine out with remaining thrust non-vectored. 



(b) Right engine out with remaining thrust vectored through eg 
Figure 4. - A possible lateral thrust vectoring method;.. 
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(a) Close view of digital equipment. 
Figure 5. - Layout of experimental facilities, 
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(b) Close view of analog equipment 
Figure 5, - Concluded. 




Figure 6. - Simulation flow chart. 
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Figure 6. - Concluded. 
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Figure 7 . - Condition 1 longitudinal variables for 
no thrust vectoring. 
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Figure 7 . Concluded 


Figure 8. Condition 1 lateral variables for 
no thrust vectoring. 
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Figure iQ. - Condition 2 lateral variables for 
no thrust vectoring. 



















Figure 10. - Concluded. 
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Figure 11. - Concluded 
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Figure 13l - Condition 4 longitudinal variables for ! 

no thrust vectoring. | 
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Figure 14. - Condition 4 lateral variables for 
no thrust vectoring. 
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Figure 15. - Condition 3 longitudinal variables for 

2 sec reaction delay/ 1 sec actuation interval 
without thrust command double. 







Figure 15. - Concluded 
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Figure 17. - Condition 3 longitudinal variables for 

2 sec reaction delay/ 1 sec actuation interval 
with thrust cotnmand double. 
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Figure 18 . 


- Condition 3 lateral variables for 

2 sec reaction delay/ 1 sec actuation interval 
with thrust command double. 
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Figure I 9 . - Condition 3 longitudinal variables for 

4 sec reaction delay/2 sec actuation interval 
without thrust command double. 
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Figure 20. ~ Concluded. 




Figure 21. - Condition 3 longitudinal variables for 

4 sec reaction delay/ 2 sec actuation interval 
with thrust command double. 

















Figure 22. -- Condition 3 lateral variables for 

4 sec reaction delay/ 2 sec actuation interval 
with thrust command double. 
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